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Anisotropic finite-size scafing of an elastic string at the depinning 
threshold in a random-periodic medium 

S. Bustingorry,^*A. B. KoltorP 

We numerically study the geometry of a driven elastic string at its sample-dependent 
depinning threshold in random-periodic media. We find that the anisotropic finite-size 
scaling of the average square width and of its associated probability distribution are 
both controlled by the ratio k = M/L'^'^'p, where i^dep is the random-manifold depinning 
roughness exponent, L is the longitudinal size of the string and M the transverse peri- 
odicity of the random medium. The rescaled average square width w'^ / L^'"^"^ displays a 
non-trivial single minimum for a finite value of k. We show that the initial decrease for 
small k reflects the crossover at ~ 1 from the random-periodic to the random-manifold 
roughness. The increase for very large k implies that the increasingly rare critical configu- 
rations, accompanying the crossover to Gumbel critical-force statistics, display anomalous 
roughness properties: a transverse-periodicity scaling in spite that <C M, and sublead- 
ing corrections to the standard random-manifold longitudinal-size scaling. Our results are 
relevant to understanding the dimensional crossover from interface to particle depinning. 



I. Introduction 

The study of the static and dynamic properties of 
d-dimensional elastic interfaces in d-|-l-dimensional 
randoni media is of interest in a wide range of phys- 
ical systems. Some concrete experimental exam- 
ples are magnetic [l}[4] or ferroelectric 5][6 domain 



walls, contact lines of liquids [t], fluid invasion in 
porous media j8[[9], and fractures 10 Tl]. In all 
these systems, the basic physics is controlled by the 
competition between quenched disorder (induced 
by the presence of impurities in the host materials) 
which promotes the wandering of the elastic object, 
against the elastic forces which tend to make the 
elastic object flat. One of the most dramatic and 
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worth understanding manifestations of this compe- 
tition is the response of these systems to an external 
drive. 

The mean square width or roughness of the inter- 
face is one of the most basic quantities in the study 
of pinned interfaces. In the absence of an external 
drive, the ground state of the system is disordered 
but well characterized by a self-affine rough geome- 
try with a diverging typical width w ^ L'^"'^ , where 
L is the linear size of the elastic object and ^eq 
is the equilibrium roughness exponent. When the 
external force is increased from zero, the ground 
state becomes unstable and the interface is locked 
in metastable states. To overcome the barriers sep- 
arating them and reach a finite steady-state veloc- 
ity V it is necessary to exceed a finite critical force, 
above which barriers disappear and no metastable 
states exist. For directed d-dimensional elastic in- 
terfaces with convex elastic energies in a D = d -I- 1 
dimensional space with disorder, the critical point 
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is unique, characterized by the critical force F = 
and its associated critical configuration fT2l. This 
critical configuration is also rough and self-affine 
such that w ^ L^dop -^jth C^gp the depinning rough- 
ness exponent. When approaching the threshold 
from above, the steady-state average velocity van- 
ishes like V ^ {F — Fc)^ and the correlation length 
characterizing the cooperative avalanche-like mo- 
tion diverges as ^ ~ (F — Fc)~'^ for F > Fc, where 
P is the velocity exponent and v is the depinning 
correlation length exponent [Ts} 16 . At finite tem- 
perature and for F <^ Fc, the system presents an 
ultra-slow steady-state creep motion with universal 



features 17 18 directly correlated with its multi- 
affine geometry 19|20 . At very small temperatures 
the absence of a divergent correlation length below 
Fc shows that depinning must be regarded as a non- 



standard phase transition 20 21 while exactly at 
F — Fc, the transition is smeared-out with the ve- 
locity vanishing as ?; ~ T'^, with i/j, the so-called 



thermal rounding exponent 22 ■ 27 



During the last years, numerical simulations have 
played an important role to understand the physics 
behind the depinning transition thanks to the de- 
velopment of powerful exact algorithms. In partic- 
ular, the development of an exact algorithm able to 
target efficiently the critical configuration and crit- 



ical force for a given sample 28 29] has allowed to ity 38 



study, precisely, the self-afBne rough geometry at 
depinning [7|[28]-31 , the sample-to-sample critical 
force distribution 32 , the critical exponents of the 
depinning transition 26 27| 33 , the renormalized 



disorder correlator 34 , and the avalanche-size dis- 
tribution in quasistatic motion [35| . Moreover, the 
same algorithm has allowed to study, precisely, the 



transient universal dynamics at depinning 36 37 



and an extension of it has allowed to study low- 
temperature creep dynamics 



20 21 



In practice, the algorithm for targeting the crit- 
ical configuration [28j[29] has been numerically 
applied to directed interfaces of linear size L dis- 
placing in a disordered potential of transverse di- 
mension M, applying periodic boundary conditions 
in both directions in order to avoid border effects. 
This is thus equivalent to an elastic string displac- 
ing in a disordered cylinder. The aspect ratio be- 
tween longitudinal L and transverse M periodic- 
ities must be carefully chosen, in order to have 
the desired thermodynamic limit corresponding to 
a given experimental realization. In Ref. i32| it 



was indeed shown that the critical force distribu- 
tion P{Fc) displays three regimes associated with 
M: (i) At very small M compared with the typical 
width L'^'^'^p of the interface, the interface wraps the 
computational box several times in the transverse 
direction, as shown schematically in Fig.[l|b), and 
therefore the periodicity of the random medium is 
relevant and P{Fc) is Gaussian; (ii) At very large 
M compared with L^<i=p , as shown schematically in 
Fig.[T]^c), periodicity effects are absent but then the 
critical force, being the maximum among many in- 
dependent sub-critical forces, obeys extreme value 
statistics and P{Fc) becomes a Gumbel distribu- 
tion; (iii) In the intermediate regime, where M « 
2^Cdop gj^fj^ periodicity effects are still irrelevant, as 
shown schematically in Fig. [l|^a) , the distribution 
function is in between the Gaussian and the Gum- 
bel distribution. It has been argued that only the 
last case, where M sa L'-'^^p, corresponds to the 
random-manifold depinning universality class (pe- 
riodicity effects absent) with a finite critical force 
in the thermodynamic limit L, Af — > oo. This cri- 
terion does not give, however, the optimal value 
of the proportionality factor between M and L^'J^p , 
and must be modified at finite velocity since the 
crossover to the random-periodic universality class 
at large length-scales depends also on the veloc- 
To avoid this problem, it has been there- 



fore proposed to define the critical scaling in the 
fixed center of mass ensemble 39 . The crossover 



from the random-manifold to the random-periodic 
universality class is, however, physically interest- 
ing, as it can occur in periodic elastic systems such 
as elastic chains. Remarkably, although the map- 
ping from a periodic elastic system (with given lat- 
tice parameter) in a random potential to a non- 
periodic elastic system (such as an interface) in 
a random potential with periodic boundary con- 
ditions is not exact, it was recently shown that the 
lattice parameter does play the role of M for elastic 
interfaces with regard to the geometrical or rough- 
Since the periodicity can often 



38 



ness properties 
be experimentally tuned in such periodic systems 
it is thus worth studying in detail the geometry 
of critical interfaces of size L as a function of M 
with periodic boundary conditions, and thus com- 
plement the study of the critical force in such sys- 
tems [32] . 

In this paper, we study in detail, using numerical 
simulations, the geometrical properties of the one- 
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Figure 1: (a) Elastic string driven by a force F in 
a random-periodic medium with periodic boundary 
conditions. It is described by a displacement field 
u{z) and has a mean width w. The anisotropic 
finite-size scaling of width fluctuations are con- 
trolled by the aspect-ratio parameter k = M/L'^'^^p , 
with Cdcp the random-manifold roughness exponent 
at depinning. In the case fc <C 1 (b) periodicity ef- 
fects are important, while when A; 3> 1 (c) they 
are not important but the roughness scaling of the 
critical configuration is anomalous. 



dimensional interface or elastic string critical con- 
figuration in a random-periodic pinning potential 
as a function of the aspect ratio parameter k, con- 
veniently defined as fc = M/L'»'^'=p. We show that fc 
is indeed the only parameter controlling the finite- 
size scaling (i.e. the dependence of observables with 
the dimensions L and M) of the average square 
width and its sample-to-sample probability distri- 
bution. The scaled average square width w^L~^'''*"p 
is described by a universal function of fc displaying 
a non-trivial single minimum at a finite value of 
fc. We show that while for small fc this reflects the 
crossover at fc ~ 1 from the random-periodic to 
the random-manifold depinning universality class, 
for large fc it implies that in the regime where the 
depinning threshold is controlled by extreme value 
(Gumbel) statistics, critical configurations also be- 
come rougher, and display an anomalous roughness 
scahng. 



II. Method 

The model we consider here is an elastic string in 
(1-1-1) dimensions described by a single valued func- 
tion u{z,t), which gives the transverse displace- 
ment u as a function of the longitudinal direction z 
and the time t [see Fig.[lja)] . The zero-temperature 
dynamics of the model is given by 

jdtu{z,t) =cd^Mz,t) + Fp{u,z)+F, (1) 

where 7 is the friction coefficient and c the elas- 
tic constant. The first term in the right hand side 
derives from an harmonic elastic energy. The ef- 
fects of a random-bond type disorder is given by the 
pinning force Fp(u,z) = —duU{u,z). The disorder 
potential U(u,z) has zero average and sample-to- 
sample fluctuations given by 



[C/(u, z) - U[u', z')Y = 5{z - z') R^{u - u'), (2) 

where the overline indicates average over disorder 
realizations and R{u) stands for a correlator of fi- 
nite range rf |18|. Finally, F represents the uni- 
form external drive acting on the string. Physi- 
cally, this model can phenomenologically describe, 
for instance, a magnetic domain wall in a thin film 
ferromagnetic material with weak and randomly lo- 
cated imperfections [Ij , being F proportional to an 
applied external magnetic field pushing the wall in 
the energetically favorable direction. 

In order to numerically solve Eq. ([T]), the system 
is discretized in the z-direction in L segments of 
size (5z = 1, i.e. z — > j = 0, L — 1, while keeping 
Uj{t) as a continuous variable. To model the con- 
tinuous random potential, a cubic spline is used, 
which passes through M regularly spaced uncor- 
related Gaussian number points [30|. For the nu- 
merical simulations performed here we have used, 
without loss of generality, 7 = 1, c = 1 and = 1 
and a disorder intensity i?(0) = 1. In both spatial 
dimensions we have used periodic boundary condi- 
tions, thus defining a L x M system. 

The critical configuration Uc{z) and force Fc are 
defined from the pinned (zero-velocity) configura- 
tion with the largest driving force F in the long 
time limit dynamics. They are thus the real solu- 
tions of 



cdlu{z) + Fp{u,z) + F ^ 0, 



(3) 
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such that for F > Fc there are no further real so- 
lutions (pinned configurations). Middleton theo- 
rems |12j assure that for Eqs. ([3| the solution ex- 
ists and it is unique for both Uc{z) and Fc, and that 
above F^. the string trajectory in an L dimensional 
phase-space is trapped into a periodic attractor (for 
a system with periodic boundary conditions as the 
one we consider). In other words, the critical con- 
figuration is the marginal fixed point solution or 
critical state of the dynamics, being Fc the criti- 
cal point control parameter of a Hopf bifurcation. 
Solving the L-dimensional system of Eqs. ^ for 
large L directly is a formidable task, due to the 
non-linearity of the pinning force Fp. On the other 
hand, solving the long-time dynamics at different 
driving forces F to localize Fc and Uc is very in- 
efficient due to the critical slowing down. Fortu- 
nately, Middleton theorems, and in particular the 
"non-passing rule", can be used again to devise a 
precise and very efficient algorithm which allows to 
obtain the critical force Fc and the critical configu- 
ration u'^ for each independent disorder realization 
iteratively without solving the actual dynamics nor 
directly inverting the system of Eqs. ^ 30j. Once 
the critical force and the critical configuration are 
determined with this algorithm, we can compute 
the different observables. In particular, the square 
width or roughness of the string at the critical point 
for a given disorder realization is defined as 



w 



3=0 



L-l 



k=0 



(4) 



Computing w'^ for different disorder realizations al- 
lows us to compute its disorder average and the 
sample-to-sample probability distribution P{w^). 
In addition, the average structure factor associated 
to the critical configuration is 



L-l 
3=0 



(5) 



where q = 2i:n/L, with n — \,...,L — 1. One 
can show, using a simple dimensional analysis, that 
given a roughness exponent such that ^ L^'», 



the structure factor behaves as S{q) 



-(1+20 

for small q, thus yielding an estimate to C with- 
out changing L. To compute averages over disor- 
der and sample-to-sample fluctuations, we consider 



i 10 




Figure 2: The scaling of w'^ for the critical con- 
figuration at different M values as indicated. The 
curves for M — 64 and 16384 are shifted upwards 
for clarity. The dashed and dotted lines are guides 
to the eye showing the expected slopes correspond- 
ing to the different roughness exponents. 



between 10'^ and 10^ independent disorder realiza- 
tions depending on the size of the system. 



III. Results 

i. Roughness at the critical point 

Figure [2] shows the scaling of the square width of 
the critical configuration with the longitudinal 
size of the system L for L = 32, 64, 128, 256, 512 and 
different values of M. When M is small, M = 8, for 
all the L values shown we observe w"^ ^ L^'^^ with 
Cl — 1.5, corresponding to the Larkin exponent in 
(1 -I- 1) dimensions. This value is different from the 
value Cdcp = 1.25 [33 40 expected for the random- 



manifold universality class, and is thus indicating 
that the periodicity effects are important for this 
joint values of M and L. This situation is schemat- 
ically represented in Fig. [l](b). This result is a 
numerical confirmation of the two-loop functional 
renormalization group result of Ref. 16 which 
shows that the C = fixed point, leading to a uni- 
versal logarithmic growth of displacements at equi- 
librium is unstable. The fluctuations are governed, 
instead, by a coarse-grained generated random- 
force as in the Larkin model, yielding a roughness 
exponent Cl = (4 — d)/2 in d dimensions 16 , which 
agrees with our result for c? = 1. We can thus say 
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Figure 3: Structure factor of the critical configu- 
ration for L = 256 and different M values, as indi- 
cated. The curves for M = 64 and 16384 are shifted 
upwards for clarity. The dashed and dotted lines 
are guides to the eye showing the expected slopes 
corresponding to the different roughness exponents. 



Figure 4: Scaling of the structure factor of the crit- 
ical configuration for L = 256 and different values 
of the transverse size M = 2^ with p = 3, 4, 14 
M . Although the values of the two exponents are 
very close, the change in the slope of the scaling 
function against the scaling variable x — q M^/^''"=p 
is clearly observed. 



that for small enough M (compared to L) the sys- 
tem belongs to the same random-periodic depin- 
ning universality class as charge density wave sys- 
tems 14 41 , which strictly correspond to Af = 1. 

When M is large, on the other hand, M — 
16384 in Fig. [2j for all the L values consid- 
ered the exponent is consistent with Cdcp, of the 
random-manifold universality class. This situation 
is schematically represented in Fig. [l|c), and we 
will show later that, for this elongated samples, the 
effects of extreme value statistics are already visi- 
ble. 

For intermediate values of A/, such as M = 
64 in Fig. [2j we can observe the crossover in 
the scale-dependent roughness exponent C,{L) ^ 

\ 'd'kfgl changing from ^dcp to Cl as L increases, 
as indicated by the dashed and dotted lines. 
This crossover, from the random-manifold to the 
random-periodic depinning geometry, occurs at a 
characteristic distance I* ~ M^/^'''=p, when the 
width in the random-manifold regime reaches the 
transverse dimension or periodicity M. At finite 
velocity, this crossover length remains constant up 
to a non-trivial characteristic velocity and then de- 



creases with increasing velocity 38 



The above mentioned geometrical crossover can 
be studied in more details through the analysis of 
the structure factor S{q), for a line of fixed size 
L. In Fig. |3] we show S{q) for L = 256 and M = 



8, 64, 16384. For the intermediate value M = 64 a 
crossover between the two regimes is visible, and 
can be described by 



^-(1+2Cl) 

g^(l + 2Cdcp) 



(6) 



with q* expected to scale as q* ^ l*^^ ^ M^^^'^'^'p. 
Therefore, the structure factor should scale as 
S'gM~'^^+^/'»'i'=p^ = H{x), where the scaled variable 
is X = gM^/^'i^p ~ q/q* and the scaling function 
behaves as 



H{x) 



2,-(1+2Cl) x^l, 

a;-(i+2Cdop) a; > 1. 



(7) 



The collapse of Fig.|4]for L — 256 and different val- 
ues of M = 2P with p = 3, 4, 14 shows that this 
scaling form is a very good approximation. How- 
ever, as we show below, small corrections can be 
expected fully in the random-manifold regime in 
the large ML^^'^^i' limit of very elongated samples. 

In Fig. [sja), we show w'^ as a function of the 
transverse periodicity M for different values of the 
longitudinal periodicity L. Remarkably, w"^ is a 
non-monotonic function of M . For small M it de- 
creases towards an L dependent minimum m* , and 
then increases with increasing M, in the regime 
where the extreme value statistics starts to affect 
the distribution of the critical force 32 . Since 
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Figure 5: (a) Squared width of the critical con- 
figuration as a function of M for different system 
sizes L as indicated, (b) Scahng of the width in 
(a), showing that the relevant control parameter 
is M/L^'''=p. The dashed line in (a) and (b) corre- 
sponds to w"^ — Af^, which is always to the left of 
the minimum of w"^ occurring at k* = m* L~'^'^^p . 
The solid line indicates k'^('^^'^i-/'^<iap) which is the 
behavior expected purely from the random-periodic 
to random-manifold crossover at the characteristic 
distance I* - MVCdop^ 



the only typical transverse scale in Fig. [sj^a) is 
set by the minimum m* , we can expect ~ 
m*^G{M/ni*) with G{x) some universal function. 
On the other hand, since the only relevant char- 
acteristic length-scale of the problem is set by the 
crossover between the random-periodic regime and 
the random-manifold regime, we can simply write 
771* ijCdop g^j^j therefore 



^2 2,-2Cdcp _ GiML-'^"'"') 



(8) 



This scaling form is confirmed in Fig. [5]jb) 
and shows that the aspect-ratio parameter k = 
A/L^^'^^p fully controls the anisotropic finite-size 
scaling of the problem. It is worth, however, not- 
ing some interesting consequences of the result of 
Fig. [sjb), as we describe below. 

Since at very small k the interface is in the 
random-periodic regime, Eq. ^ should led to w"^ ~ 
L'^'^^ and therefore one deduces that. 



G(fc) - fc^^i-^i'/'^'i-p), fc<fc* 



(9) 



where k* — m*L^^'^<=p. The fact that the random- 
periodic roughness exponent Cl = 3/2 is larger than 
the random-manifold one ^dep ~ 5/4 consequently 
implies an initial decrease of G(fc) as G{k) ~ 
as shown in Fig. [5jb) by the solid line. Periodicity 
effects, or the crossover from random-periodic to 
random-manifold, thus explain the initial decrease 
of G{k) seen in Fig. [sjjb), or the initial decrease of 

against M for fixed L, seen in Fig. (sj^a). At 
this respect, it is then worth noting that the line 
w'^ = M^, shown by a dashed line, lies completely 
in the regime k < k* implying that the naive cri- 
terion < is not enough to avoid period- 
icity effects, and to have the system fully in the 
random-manifold regime. As we show later, this 
is related with the shape of the probability distri- 
bution of P{w^) which displays sample-to-sample 
fluctuations of the order of the average w^. 

The presence of a minimum at k* in the function 
G{k) and in particular its slower than power-law 
increase for A; > fc* is non-trivial and constitutes 
one of the main results of the present work. This 
result shows that corrections to the standard scal- 
ing w'^ ~ iCdep jj^oy arise from the aspect-ratio de- 
pendence of the prefactor G{k). On the one hand, 

now grows with M for L fixed, in spite that 
w"^ <C M^, i.e. transverse-size/periodicity scaling is 
present. On the other hand, the scaling of with 
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L is slower in this regime, due to subleading scaling 
corrections coming from G{k). The precise origin 
of these interesting leading and subleading correc- 
tions in the finite-size anisotropic scaling are highly 
non-trivial. Since the critical configurations in this 
regime have the constant roughness exponent Cdep 
of the random-manifold universality class, the slow 
increase of G{k) cannot be attributed to a geomet- 
rical crossover effect, as for the case k < k* . How- 
ever, we might relate this effect to the crossover 
in the critical force statistics, from Gaussian to 
Gumbel, in the k :$> k* limit [32]. In the Gum- 
bel regime, the average critical force is expected to 



increase as Fr 



39 



smce 



log(M/L^3p) = logfc _ 
the sample critical force can be roughly regarded 
as the maximum among Af/i^ep independent sub- 
critical forces and configurations (32l . The increase 
in the critical force might be therefore correlated 
with the slow increase of roughness. The physi- 
cal connection between the two is subtle though, 
since a large critical force in a very elongated sam- 
ple could be achieved both by profiting very rare 
correlated pinning forces such as accidental colum- 
nar defects, or by profiting very rare non-correlated 
strong pinning forces. Since in the first case the 
critical configuration would be more correlated and 
in general less rough than for less elongated sam- 
ples (smaller k), contrary to our numerical data of 
Fig. [sjjb) , we think that the second cause is more 
plausible. We can thus think that in the k ^ k* 
limit of extreme value statistics of Fc, the effective 
disorder strength on the critical configuration in- 
creases with k. This might be translated into the 
universal function G{k), such that w-^ « L^'''^'=pG{k) 
can increase for increasing values of k at fixed L in 
such regime. A quantitative description of these 
scaling corrections remains an open challenge. 

ii. Distribution function 

We now analyze sample-to-sample fluctuations of 
the square width w'^ by computing its probability 
distribution P{w'^). This property is relevant as 
fluctuates even in the thermodynamic limit for 
critical interfaces with a positive roughness expo- 
It has been computed for models with 

or 



nent 



42 



43 



44 45 , the Mulhns 



dynamical disorder such as random-walk 
Edwards- Wilkinson interfaces 
Herrings model [46] and for non-Markovian Gaus- 
sian signals in general 47 48 . It has also been 



1.2 
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Figure 6: Scaling function $(a;) for L — 256 and 
different values of M 8, 128, 2048, 16384, which 
shows the change with the transverse size M. 



calculated for non-linear models such as the one- 
dimensional Kardar-Parisi-Zhang model 



49 50 



and for the quenched Edwards- Wilkinson model at 



equilibrium 51 



In particular, the probability distribution P{w^) 
of critical interfaces at the depinning transition 



was studied analytically 52 , numerically 31 and 



also experimentally for contact lines in partial 
wetting [?]. Remarkably, non-Gaussian effects in 
depinning models are found to be smaller than 
0.1% [3l]|52], thus showing that P{w'^) is strongly 
determined by the self-afhne (critical) geometry it- 
self, rather than by the particular mechanism pro- 
ducing it. As in all the above mentioned systems 
the width distribution P(w^) at different universal- 
ity classes of the depinning transition was found to 
scale as 



(10) 



with an universal function, which only depends 
on the roughness exponent C and on boundary con- 
ditions when the global width is considered iTfiS 



In this way, w"^ is the only characteristic length- 
scale of the system, absorbing the system longitu- 
dinal size L, and all the non-universal parameters 
of the model such as the elastic constant of the 
interface, the strength of the disorder and/or the 
temperature. Since $^ can be easily generated us- 
ing non-Markovian Gaussian signals '53] , the quan- 
tity w'^PiuP') is a good observable to extract the 
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Figure 7: Scaling function for different val- 

ues of L = 32,64,128,256 while keeping (a) k = 
M/L<^-^^ w 1 and (b) k = M/L'^^-" w 0.025. 
The dotted line corresponds to the scaling function 
of the non-disordered Edwards-Wilkinson equa- 
while the continuous and dashed lines cor- 



43 



tion 

respond to the scaling functions of Gaussian signals 
with ^ = 1.25 and ( — 1.5, respectively (31„53, . 



Fig. ^a) and k « OMb < k* in Fig. [Tjjb), with 
k* the minimum of w^. Since data for the same 
k practically collapses into the same curve, we can 
write for our case: 



w 



(11) 



Therefore, the anisotropic scaling of the probability 
distribution is fully controlled hy k, as it was found 
for w^. 

In Figs, and (b), we also show the uni- 

versal functions and ^^^^p generated using 

non-Markovian Gaussian signals 31 53 , and for 
comparison we also show $1/2 corresponding to 
the Markovian periodic Gaussian signal or the 
Edwards- Wilkinson equation 43 . Comparing this 
with the collapsed data for depinning, we see that 

the function <i> ( SL A; ) respects the limits 



$(a;,/c^O) 
> k*) 



(12) 



as expected from the existence of the geometric 
crossover between the roughness exponents Cl for 
fc — > and Cdcp for k > k* . For intermediate values 

k < k* , however, does not necessarily 

\ w-^ J 

coincide with the one of a Gaussian signal func- 
tion $^ for a given C, since the critical configura- 
tion includes a crossover length 1* < L. Whether 
multi-affine or effective exponent self-affine non- 
Markovian Gaussian signals can be used to describe 
satisfactorily these intermediate cases is an inter- 
esting open issue. 



roughness exponent of a critical interface from ex- 
perimental data. 

In Fig. § we show how the scaled distribution 
function ^(x) = iiP- P{x ■ufl') looks like for the de- 
pinning transition in a random-periodic medium for 
a fixed value L = 256 and different values of M. We 
see that $(x) depends on M for small M but con- 
verges to a fixed shape for large M. We also note 
that for all M $(a;) extends appreciably beyond 
X = \ explaining why the criterion uP- < is not 
enough to be fully in the random-manifold regime, 
as noted in Fig. [5j 

In Fig. [rj we show the scaling function $(x) for 
different values of L and M but fixing the aspect- 
ratio parameter k = M/L^'^'^'p, fc « 1 > fc* in 



IV. Conclusions 

We have numerically studied the anisotropic finite- 
size scaling of the roughness of a driven elastic 
string at its sample-dependent depinning thresh- 
old in a random medium with periodic bound- 
ary conditions in both the longitudinal and trans- 
verse directions. The average square width w'^ and 
its probability distribution are both controlled by 
the parameter k = M/L'^'^'^^. A non-trivial sin- 
gle minimum for a finite value of k was found in 
u'2/L^f''''p . For small fc, the initial decrease of w"^ re- 
flects the crossover from the random-periodic to the 
random-manifold roughness. For very large fc, the 
growth with k implies that the crossover to Gumbel 
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statistics in the critical forces induces corrections 
to G{k), that grow with k, to the string rough- 
ness scaling w"^ ~ G{k)L'^'^'^<'^ . These increasingly 
rare critical configurations thus have an anoma- 
lous roughness scaling: they have a transverse- 
size/periodicity scaling in spite that its width is 
w"^ <C M^, and subleading (negative) corrections 
to the standard random-manifold longitudinal-size 
scahng. 

Our results could be useful for understanding 
roughness fluctuations and scaling in finite ex- 
perimental systems. The crossover from random- 
periodic to random-manifold roughness could be 
studied in periodic elastic systems with variable 
periodicity, such as confined vortex rows ,54 and 



single-files of macroscopically charged particles 55 
or colloids 



56 



with additional quenched disor- 
der. The rare-event dominated scaling corrections 
to the interface roughness scaling could be studied 
in systems with a large transverse dimension, such 



as domain walls in ferromagnetic nanowires 57 



For the later case, it would be interesting to have 
a quantitative theory, making the connection be- 
tween the extreme value statistics of the depinning 
threshold and the anomalous scaling corrections to 
the roughness of such rare critical configurations. 
This would allow to understand the dimensional 
crossover, from interface to particle depinning. 
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